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Our result in Ref. yjj that, instead of the con- 
ventional heat fluctuation theorem (FT), a new FT 
holds for heat fluctuations for a Brownian particle 
in a moving confining potential^, was claimed to be 
disproved in a very recent preprint by Narayan and 
DharWI/. This comment is meant to show that their 
assertion is not correct. The point is that they formu- 
late their FT differently than we do ours. Effectively, 
their FT speaks about a physically irrelevant limiting 
case of our new FT. This implies that the two FTs 
are not in contradiction with each other. Further- 
more, we point out an incorrect assumption in their 
derivation. 

First, let us remark that although we have a FT 
for all times in Ref. pQ, we only discuss the infinite 
time version in this comment. 

What we consider as the conventional and the new 
FT, and what Narayan and Dhar in Ref. [3] consider 
the FT to be, respectively, can be formulated as fol- 
lows. 

We pQ: P(Q T — pcrr) denotes the probability that 
the heat Q T produced in a time interval r, has the 
value pot, where a is the average heat produced per 
unit time (= (Q t )/t) in the stationary state over 
positive times0J|5].* 

Narayan and Dhar 3 : P(Q T — A) denotes the 
probability that the heat Q T produced in a time r 
has the fixed value A, where A is independent of r. 
We note that in Ref. |3j this is formulated in terms of 
S, which is the entropy production in a time r, and 
is given by S = Q T /T, with T the temperature of the 
fluid in which the Brownian particle moves. 

As a consequence, there are two similarly look- 
ing FTs, but at closer inspection they are physically 



radically different. 

We PP consider the conventional FT in the form 
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for fixed p, where p characterizes Q T in units of 
err = (Q T ) and (3 = 1/kgT. This is consistent with 
Refs. Ill EJ. 

Narayan and Dhar [2] consider the conventional 
FT to be 
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for fixed, finite A. 

Narayan and Dhar's FT is so formulated — be- 
cause of the independence of A from r — that when 
r approaches infinity, so that (Q T ) also approaches 
infinity, the fluctuation value A remains always fi- 
nite and therefore much smaller than (Q T ). Then 
only values of A much smaller than (Q T ), but no val- 
ues either near or larger than (Q T ), are considered in 
Narayan and Dhar's FT. Hence, only an increasingly 
restricted range of fluctuation values A are actually 
studied in this FT as r approaches infinity. 

We, in contrast, consider all fluctuation values A, 
because we scale the value A with r, so that when 
(Q T ) grows (~ r), so does A — par, as is written 
explicitly in our formulation in Eq. for the FT. 
A r-proportional A means that for a given (Q T ), any 
heat fluctuations (Q T — par) can be considered by 
an appropriate choice of p, even as r approaches in- 
finity J 

So we claim that Eq. (JTJ is consistent with Eq. © 
since a constant finite A — as Narayan and Dhar's FT 



*In Ref. d, cr is called w. 

TThis issue of scaling is analogous to using the proper scaling to derive, e.g., the diffusion equation for the random walk, 
which requires an I 2 /t scaling. Other scalings (like e.g., l/t) would not lead to a physically meaningful description of diffusion 7 . 
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uses in Eq. — can a ^ so be obtained from Eq. QJ, 
by taking the limit p — > (when r — ► oo). In other 
words, Eq. J2J is a very special case of Eq. which 
covers strictly speaking only the value p = 0, while 
Eq. covers all values of p. 

We will now argue that the FT of Narayan and 
Dhar in Eq. (J2J is also consistent with our new FT 
found in Ref. pQ. In Ref. PP the conventional FT 
[Eq. Q] was found to be incorrect as a FT for heat 
fluctuations in the system described above, and had 
to be replaced by a new FT. However, for the special 
region —1 < p < 1, this new FT turns out to coincide 
with the conventional form of Eq. (JTJ , which includes 
p = 0. Therefore, there is no contradiction between 
our new FT and Narayan and Dhar's FT in Eq. © 
either. 

The above is enough to conclude that our new 
FT in Ref. pQ is not in contradiction with the result 
Eq. J5J (or their Eq. (1) in Ref. [3])- Moreover, our 
theory leading to the new FT in Ref. pQ has been 
compared in detail with two independent numerical 
calculations, one, a numerically carried out inverse 
Fourier transform explained in Ref. PP, and another, 
a sampling method, which is yet to be published jo]. 
Both give perfect agreement with each other as well 
as with the new FT. So also from this, we have no 
reason to distrust the results in Ref. pQ . 

We take this opportunity to address very briefly 
two other points in Ref. 

Firstly, in their derivation of their (conventional) 
FT for Q T in their Sec. Ill, Narayan and Dhar ne- 
glect a term in their Eq. (17) which one can neglect 
for the average (Q T ), but not for the fluctuations of 
Q T . In fact, it is precisely this term which gives rise 
to exponential rather than Gaussian behavior of the 
tails of the distribution P{Q T ) — the existence of 
which Narayan and Dhar also notice — and through 



these exponential tails, this term also gives rise to our 
new FT, instead of to the conventional one. To sub- 
stantiate this, a saddle point method|TJ|H|, involving 
a complicated study of the behavior of the Fourier 
transform of P(Q T ), in particular, its singularities in 
the complex plane, is needed, which goes beyond the 
scope of this comment. 

Secondly, Narayan and Dhar seem to suggest that 
the singularity of the Langevin equation for vanish- 
ing mass m of the particle might be a problem in the 
derivation of a FT. However, setting m = from the 
start in their Eq. (12) or setting m = at the end 
in their Eqs. (17) and (19), give the same FT, i.e., 
no singularity is encountered in this derivation with 
respect to m = 0. 
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